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Abstract 

We explore the parameter space of a Constrained Next-to- Minimal Supersymmetric 
Standard Model with GUT scale boundary conditions (CNMSSM), and find regions 
where the relic density of the lightest neutralino is compatible with the WMAP mea- 
surement. We emphasize differences with the MSSM: cases where annihilation of the 
LSP occurs via a Higgs resonance at low values of tan/3 and cases where the LSP has 
a large singlino component. The particle spectrum as well as theoretical and collider 
constraints are calculated with NMSSMTools. All neutralino annihilation and coannihi- 
lation processes are then computed with micrOMEGAs, taking into account higher order 
corrections to the Higgs sector. 

1 Introduction 

One of the attractive features of supersymmetric (susy) extensions of the Standard Model 
(SM) with conserved R-parity, is the presence of a good Dark Matter (DM) candidate, the 
lightest SUSY particle (LSP). Over the years, numerous studies have examined the con- 
straints on the parameter space of the constrained minimal supersymmetric standard model 
(CMSSM) including that originating from the cold DM abundance extracted from preci- 
sion cosmological measurements, notably those of WMAP [1] and SDSS [2]. All concluded 
that the CMSSM could in some region of parameter space provide a satisfactory DM can- 
didate [3-10]. The MSSM, or its constrained version, however face a naturalness problem 
- the so-called /i problem. The NMSSM is a simple extension of the MSSM that solves 
this problem elegantly via the introduction of a gauge singlet superfield, S. The effective 
yU parameter term is determined by the vev of this singlet field, which is naturally of EW 
scale [11-19]. 

The NMSSM contains an extra scalar and pseudoscalar states in the Higgs sector as well 
as an additional neutralino, the singlino. Owing both to modifications in the Higgs sector 
and the neutralino sector, the DM properties can differ from those of the MSSM [20-25]. In 
particular a LSP with a large singlino component has different annihilation properties than 
the bino LSP that is in general found in the CMSSM. It is even possible to have a very 
light singlino LSP if accompanied by a very light scalar. Even when the LSP has no singlino 
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component, the more elaborate Higgs sector of the model provides additional channels for 
rapid annihilation through Higgs exchange. This can have implications for direct and indirect 
detection rates [25-28]. 

The purpose of this paper is to explore a constrained version of the NMSSM (CNMSSM) 
with semi-universal parameters defined at the GUT scale. Our choice of semi- universality 
(i.e. non universal singlet sector) rather than strict universality is motivated in sec. 2. In 
this first analysis we emphasize the differences with the CMSSM predictions as concerns 
the neutralino properties and its annihilation. We will therefore concentrate on regions of 
parameter space where the LSP has some singlino component as well as regions where the 
spectrum is such that annihilation can take place near a Higgs resonance. In the former case, 
we will see that coannihilation processes, with staus or other neutralinos play a crucial role. 
Of course, in large regions of parameter space we expect to recover features of the CMSSM 
with a LSP that is mostly bino and can only annihilate efficiently near a Higgs resonance 
or coannihilate with light sleptons. In this case the only difference with the CMSSM will 
be additional allowed parameter space due to relaxed constraints from LEP on the Higgs 
sector. 

To evaluate the supersymmetric spectrum we use the NMSPEC program from the 
NMSSMTools package [29]. Using renormalization group equations (RGBs) and starting from 
GUT scale parameters, this code computes the Higgs spectrum including higher order cor- 
rections as well as the masses of sparticles at one-loop. It also checks all the available collider 
constraints as well as various theoretical constraints. For the computation of the relic density 
of DM we rely on micrOMEGAs [30] which is included in NMSSMTools. 

The paper is organized as follows: in sec. 2 we briefly describe the model. In sec. 3 
we discuss the main channels for annihilation. In sec. 4 we present typical case studies. 
Conclusions follow in sec. 5. 



2 The CNMSSM 

In the present paper we discuss the NMSSM with a scale invariant superpotential 

W = XSHuHd +'^S^ + (Yukawa couplings) , (2.1) 

where the weak scale originates from the soft SUSY breaking scale only, i.e. where no super- 
symmetric dimensionful parameters as /i are present in the superpotential. The soft SUSY 
breaking terms in the Higgs sector are then given by 

Vsoit = mjfjHul'^ + mjfjHdl'^ + ml\S\'^ 

+ (^A^H^H^S + ^-kA^S^ + h.c.^ . (2.2) 

A possible cosmological domain wall problem [31] caused by a global .^3 symmetry can be 
avoided by introducing suitable non-renormalizable operators [32, 33] which neither generate 
dangerously large singlet tadpole diagrams [11,34-36] nor affect the low energy phenomeno- 
logy. Other models solving the MSSM /i problem with extra gauge singlet flelds and with no 
domain walls include the nearly minimal SUSY model (nMSSM), with an additional tadpole 
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term for the singlet in the superpotential and/or in the soft scalar potential [37-39], and 
the UMSSM, with an extra U{iy gauge symmetry [40-43]. The DM properties as well as 
the possibility of generating the baryon asymmetry of the universe within these models have 
been analysed in refs. [44-53]. 

Constraining the parameters of the NMSSM by imposing universality at the GUT scale 
is not as direct as in the MSSM. In the CMSSM, the free parameters are the universal GUT 
scale soft terms Aq, tuq and Mi/2- In addition, Mz and tan/5 (at the weak scale) are used as 
inputs, and the two minimization equations of the Higgs potential w.r.t. the two real Higgs 
vevs hu and are used to compute fi and B in terms of the other parameters (this leaves 
the sign of as a free parameter). Both /i and B have only a small effect on the RGEs of 
the other parameters (via threshold effects from particles whose masses depend on fi and/or 
B). In numerical codes, this is usually solved by an iterative procedure. 

At first sight, an application of this procedure to the NMSSM is not possible: neither 
fi nor B are present and one has to cope with three coupled minimization equations w.r.t. 
the Higgs vevs hu, hd and s. This means that starting from strict universality, with M1/2, 
mo, Aq as well as A and k as free parameters at the GUT scale, one usually ends up with 
the wrong value of Mz at the weak scale (as no dimensionful parameter is left to tune the 
correct value). In addition, tan/3 cannot be a free parameter in this approach. To solve this 
problem, the first phenomenological studies of the NMSSM [16-18,54] used the fact that 
Ml/2 can be factorised out of the RGEs and took universal "scaled" parameters at the GUT 
scale as inputs: Aq/Mii2 and mo/Mi/2, as well as A and k. In this approach however, tan/? 
and Ml/2 are output parameters, computed from the minimization of the potential and the 
known value of Mz- Furthermore, one usually finds that in the strictly universal NMSSM 
AQ.mQ <^ Ml/2 [54]. It is therefore difficult to extend the CMSSM DM studies, usually 
presented as plots in the mo,Mi/2 plane for fixed tan/?. In this paper we follow the semi- 
universal approach first used in the program NMSPEC [29], that we shall briefly explain 
now. 

First, let us assume that A as well as all the soft terms (except m|) are known at the 
weak scale (e.g. after integration of the RGEs down from the GUT scale). One can define 
effective (s dependent) parameters at the weak scale: 



(In the following, we will often use u rather than k). The minimization equations w.r.t. 
and hd can then be solved for the effective and S, as in the MSSM, in terms of the other 
parameters (inch Mz and tan/3). From /x and B one then deduce (for A and A\ given) both 
s and K. Finally, from the minimization equation w.r.t. s, one can easily obtain the soft 
singlet mass m| in terms of all other parameters. At tree level, the minimization equations 
giving yU (up to a sign), B and m| (i.e. k, s and m|) read: 



IJ = Xs , 



V = KS ^ 



B = Ax + iy. 



(2.3) 



171% — mi, tan^/? 1 „ 
tan2/3 - 1 2 ^ 



m% 



B 



X'^{A, + 2u) - + 2u) - \\hl + hi) . 




{mj,^+mj,^ + 2ft' + \'{hl + hl)) , 



(2.4) 
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The radiative corrections to the scalar potential show a weak dependence on s, k and 
m| which can be included in the minimization equations. These become non-linear in the 
parameters to solve for and have therefore to be solved iteratively. The derived parameters 
K and m| affect the RGEs of the other parameters not only through threshold effects around 
-^susy, but also through the (3 functions. However, the numerical impact is relatively small 
such that an iterative procedure converges quite rapidly again. 

As m| is computed from the minimization equations, it is difficult to find parameters 
such that it assumes the same value as the Higgs doublet (or other scalar) soft masses 
squared at the GUT scale. On the other hand, the mechanism for the generation of soft 
SUSY breaking terms could easily treat the singlet differently from the other non-singlet 
matter multiplets [55]. Hence, we assume the following free parameters for the CNMSSM: 

• tan/5, sign(yu) at the weak scale; 

• A at the SUSY scale; 

• universal soft terms M1/2, ruQ and Aq at the GUT scale. Exceptions are: 

— m| at the GUT scale which is an output, as described above; 

— at the GUT scale which is considered as an independent parameter. 

The reason for considering at the GUT scale as a free parameter is twofold: first, 
Af^ = Aq usually leads to negative mass squared in the Higgs sector. Second, if an underlying 
mechanism for the generation of the soft SUSY breaking terms treats the singlet differently 
from the other matter fields (as it is already assumed for this will also affect the 

coupling A^ which involves the singlet only. Hence, in our semi-universal approach, k, s and 
m| are computed from the minimization equations. This implies important differences with 
the general NMSSM where all the soft parameters as well as A, k and tan/3 are taken as free 
parameters at the weak scale, especially in the singlet neutralino (singlino) sector as we shall 
see. 

Not all choices of parameters are allowed in the CNMSSM. Some lead to negative mass 
squared for scalar fields (Higgs or sfermions), others to Landau poles below the GUT scale 
for the dimensionless couplings. The NMSSMTools package also includes all the available 
experimental constraints from LEP and Tevatron on sparticle and Higgs searches (for details 
on the exclusion channels see refs. [29,56]). It is the aim of this paper to find out which 
regions in the parameter space of the CNMSSM can fulfill all the theoretical and experimental 
tests and at the same time provide the correct amount of DM. In sec. 4, we will present 
quantitative results, but let us first have a qualitative approach here. 

In order to do this, it is helpful to have a look at the symmetric 3x3 CP even Higgs 
mass matrix in the basis {Hu, Ha, S) at tree level: 



2XKfi - Xhd{Ax + 2u) 



'U 



iT'd 



Ml 



2Xhdfi - XK{Ax + 2iy) 



(2.5) 



AM,^ + u{A^ + Au) 
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To a good approximation, the 2x2 doublet subsector is diagonalized by the angle (3 which 
gives a light eigenstate h with mass 



nih = + ^ sm^2/3^ M| (2.6) 

and a heavy eigenstate H with a mass ttlh ~ close to the MSSM-like CP odd state 
(the larger m^, the better this approximation). In the NMSSM, one can define as the 
diagonal doublet term in the CP odd 2x2 mass matrix after the Goldstone mode has been 
dropped. At tree level, it has the same expression as in the MSSM: 

It is already well known that large values of A are not allowed in the NMSSM as they 
would imply a Landau pole below the GUT scale [11-19]. This leads to an upper bound on 
A < .7 which depends on tan/3. Values of A > .1 are not always allowed in the CNMSSM: on 
the one hand A increases the mass of the light state h. This increase is relevant, however, only 
for small tan/3. On the other hand, A induces mixings between the doublet and the singlet 
states. Since the singlet state is typically heavier than ~ 100 GeV, this mixing reduces the 
mass of the lightest eigenstate h which will then often violate bounds from LEP. In order 
to maximize the mass of the lightest CP even mass eigenstate, this singlet-doublet mixing 
has to vanish (as described in [57]), which implies a relation between u, Ax, A and tan/3. 
This relation is generally not satisfied within the CNMSSM; then - at least for large values 
of tan/3 - this mixing effect disallows values of A > .1. 

When A is small, the mixings between the scalar and pseudoscalar singlet and the doublet 
states are small and their masses read, respectively: 

m| = z/(y4« + 4z/) , nip = -3uA^ . (2.8) 

The parameter being only slightly renormalized from the GUT scale down to the SUSY 
scale, eq. (2.8) shows that the masses of the singlet states are directly proportional to the 
value of at the GUT scale. The condition that both squared masses are positive together 
with eq. (2.3) implies 

-4(5 -Aa)' < A^iB-Ax) < 0. (2.9) 

The parameter B, obtained from the minimization equations (2.4) depends on mo, M1/2, A 
and tan/3, while Ax depends on Aq, A and tan/3. This means that, for sign(/i) positive (which 
we will always assume in the following), either Ai^ > and Aq > ylo(mo, M1/2, A, tan/3) or 
Af^ < and Aq < /lo(mo, M1/2, A, tan/3). Moreover, for A^ > and tan/3 moderate, large 
values of niQ or M1/2 (implying B large and positive) lead to a negative mass squared for 
the pseudoscalar singlet and are therefore disallowed. If tan/3 is large, however, B remains 
small up to large values for niQ and M1/2 which are no longer excluded. 
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Finallj^in the neutrahno sector, we have a symmetric 5x5 mass matrix, given in the 
basis {B,W,Hi,H2,S) by: 



Ml 



M2 - 



giK 

" V2 



9ihd 

g2hd 

V2 



(2.10) 







-/i -Xhd 
-Xhu 
2p 



When A <^ 1, the singhno decouples from other neutrahnos and its mass is rrig = 2u. It is 
difficult to guess the input parameters that lead to a singlino LSP as z/ is a derived parameter. 
Nevertheless, using eq. (2.3) one can rewrite m-^ = 2{B — A\). One can see from eq. (2.4) 
that for large values of tan/5 B is small and the mass of the singlino is simply rrig ~ —2 Ax, 
i.e. it depends mainly on Ao and is insensitive to mo and M1/2. Hence, the singlino will be 
the LSP for large values of M1/2 (where the bino is heavy). For moderate values of tan/?, the 
singlino mass depends also on B, which grows like mo and M1/2 and has the same sign as /i 
(assumed positive). Therefore, if A^ < and Aq < Aq, v is positive and a singlino LSP is 
likely to appear at small values of mo and Mi/2- On the contrary, if > and Aq'>^ A^, v 
is negative and the singlino will be the LSP at large mo or M1/2, where i/ — 0. (We remind 
that in this case, mo and M1/2 are bounded from above in order to have a positive mass 
squared for the pseudoscalar singlet). For A > .1 the singlino mixes with the higgsinos and 
it is more difficult to make qualitative statements. 



The relic density calculation follows the usual procedure of evaluating the thermally averaged 
cross section for annihilation and coannihilation of the LSP and solving for the density evo- 
lution equation numerically. To do so we use micrOMEGAs [30] adapted for the NMSSM [23]. 
For this a model has been specified in the CalcHEP [58] format with the help of LANHEP [59]. 
Given the set of input parameters, as specified in the SLHA2 format [60,61], the code then 
finds the LSP before generating with CalcHEP all the matrix elements of all relevant processes 
of annihilation and, when necessary, coannihilation. An automatic procedure for looking for 
s-channel poles is incorporated into the program such that a more precise integration rou- 
tine can be used in the event one is close to a pole. One issue that has to be treated with 
special care is the one of radiative corrections to the Higgs bosons masses and couplings. 
Those are calculated in NMSSMTools and can be very important. Since knowing the precise 
value of the Higgs mass is often crucial to a relic density calculation, we want to use the 
loop-corrected masses even though CalcHEP computes only tree-level matrix elements. To 
deal with this, as was described in [23], we write a general effective potential involving two 
doublets and a singlet which includes 10 effective parameters. Taking the radiatively cor- 
rected Higgs masses and mixings angles provided my NMSSMTools, we extract the value of 
these effective parameters. It is then a simple matter to derive the corresponding trilinear 
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and quartic couplings of the scalar sector. Note however that the effective potential includes 
only leading operators, and that in practice the number of parameters that we can extract 
from observables (masses and mixings) in the Higgs sector is limited. While these operators 
should include the dominant corrections to the masses and vertices, this procedure can fail 
or show too large sensitivity to a given parameter. A signal that there might be a problem 
with the procedure is that one of the dimensionless effective operators is larger than 1. A 
warning is issued in this case ^. 

A complete calculation of the neutralino DM in the NMSSM defined at the weak scale 
has shown that often the same mechanisms as in the MSSM for neutralino annihilation are 
at work. In the constrained version of the models, those can be classified as 

• Annihilation of a bino LSP through light sfermion exchange; this occurs at small values 
of mo, Ml/2. 

• Annihilation near a Higgs resonance (heavy or light Higgses); this requires some hig- 
gsino component. The light Higgs resonance is found for values of M1/2 near/below 
the LEP exclusion bound while the heavy Higgs resonance is found at large values of 
tan/3. 

• Annihilation of a mixed bino-higgsino LSP; this occurs at very large rriQ in the so-called 
focus point region. 

• Coannihilation of a bino with sfermions; this occurs at small values of mo when the 
stau is the NLSP or in a very small region at low mo, M1/2 where the stop is the NLSP. 
The latter is possible only for large negative values of Aq. 

Nevertheless in the NMSSM one finds important differences with the MSSM: the extra 
Higgses open up the possibility of more resonant annihilation and the singlino component 
of the neutralino LSP can alter the prediction for annihilation cross-sections. The new 
mechanism that can provide a DM candidate compatible with the WMAP results are 

• Annihilation near a pseudoscalar singlet resonance, occurring at any values of tan/3 
provided A > .1. 

• Coannihilation of a singlino LSP with sfermions for small values of mo and A <^ 1. 

• Coannihilation of a singlino LSP with higgsino NLSP; this occurs at large mo, A <^ 1 
and is more likely for large values of tan/3. 

• Coannihilation of a singlino LSP with bino NLSP, the bino rapidly annihilating through 
a Higgs resonance, for large tan/3 and A ^ 1. 

In the following section we will give specific examples where the old and new mechanisms 
for LSP (co) annihilation are at work. 

^It also means that there might be a mismatch in the Higgs to Higgs partial widths between NMSSMTools 
and micrOMEGAs since the set of radiative corrections to the Higgs masses that have been calculated is much 
more complete than the ones to the partial widths. In general the difference is < 10%. For the purpose of 
computing the relic density, the Higgs to Higgs decays can play a role only near a resonance where the total 
width is the relevant parameter. A correction to one partial width is therefore not so crucial. 
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4 Results 



We first performed a general scan over tlie parameter space of the CNMSSM for fixed values 
of tlie SM parameters in order to find regions satisfying all collider constraints and compatible 
with WMAP, that is nh"^ < .129 2. We quote only the upper bound since we assume that the 
neutralino does not necessarily account for all the DM. We used the central value of the top 
quark mass measured at the Tevatron, rrit = 171.4 GeV [62], and choose mf,(mb) = 4.214 GeV 
and as{Mz) = .1172. Note that the top quark mass does affect the value of the light Higgs 
mass as well as the calculation of the spectrum at large mo while affects the Higgs masses 
and couplings especially at large tan/3. As mentioned in sec. 2, we assumed sign(/i) > 0. 

We found 2 disconnected regions in the parameter space: One at small tan/5 ~ 2 and 
large A ~ .5, where the lightest Higgs boson is heavy enough to pass the LEP constraints 
due to the specific NMSSM tree level contribution to its mass as shown in eq. (2.6). The 
other allowed region is for small A < .1 and tan/5 > 4. In this case, the singlet sector is 
almost decoupled from the rest of the theory and the bound on the lightest Higgs doublet 
implies a lower bound on tan/? as in the CMSSM. 

We then picked 7 couple of values for A and tan/3: A = .5, tan/3 = 2 in the first region 
and A = .1 or .01, tan/3 = 5, 10 or 50 in the second. For each couple of values of A, tan/3 we 
scanned randomly on the remaining free parameters, namely mo, M1/2, ^o? ^k, in order to 
find regions in the parameter space allowed by all theoretical and experimental constraints. 
As mentioned in sec. 2, we found either < and Aq < 74o(A, tan/3) or > and 
Aq > Ao(A, tan/3), the latter case appearing only if A = .01 and/or tan/3 = 50 (the other 
values of A, tan/? always lead to light states in the Higgs sector excluded by LEP when 
A. > 0). 

Finally, we identified the values of Aq and A^ for which the main neutralino annihilation 
channel is a Higgs singlet resonance, or the LSP is mainly singlino. For the former case we 
found either A = .5, tan/3 = 2 and Aq ~ A^ large and negative, or A = .1, large Aq < and 
small Af^ < 0. For A = .1, tan/3 = 50 we also found singlet resonances at large Aq > and 
small > 0. The singlino LSP scenario on the other hand appears only for A = .01. It 
requires A^ < and Aq < Ao(tan/3) or A^ > and Aq > Ao(tan/3). We will now present 
plots in the mo, M1/2 plane for selected values of A, tan/3, Aq and A^. 

4.1 Large A: singlet resonances 

We first consider the cases where A = .5 or .1. Here, the LSP is mainly bino and is not 
expected to have a large singlino component. The constraints from Higgs searches at LEP 
can be very important, especially at small values of tan/3, when considering the rather low 
value of the top quark mass now measured at Fermilab, m^ = 171.4 GeV. 

^The most recent WMAP and SDSS results give tighter constraints on the relic density of DM [63,64]. 
However adopting a more conservative approach when fitting this data, that is allowing extra degrees of 
freedom in the cosmological model used [65], leads to the range .094 < ilh^ < .136 [66], giving results similar 
to the ones we discuss below. 
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m„[GeV] ^gjj ^ m„[GeV] ^ (1^^ 



Figure 1: (a) The WMAP allowed region (green) in the mo, plane for A = .5, tan/? = 2, 

Aq = —1300 GeV, = —1400 GeV. We show the region excluded by theoretical constraints or by 
LEP searches on sparticles (black), the region where a sfermion - the lightest stau here - is the LSP 
(blue), the LEP limit from Higgs searches (red/hatch) and the contour nih = 111 GeV (pink/hatch), 
(b) The DM relic density Qh'^ (black), the pseudoscalar singlet mass mp (red) and twice the bino 
LSP mass m^, as a function of niQ for the same choice of parameters and M1/2 = 600 GeV. 

4.1.1 Small tan/3 

The case tan/3 = 2 is very characteristic of the NMSSM. Indeed, when A is near its maximal 
value and tan/? is small, the light Higgs doublet can be heavier than in the MSSM, as already 
mentioned, so it is possible to satisfy the LEP constraints. Nevertheless, for tan/3 = 2 
and A = .5, large regions of parameter space are excluded by experimental or theoretical 
constraints. For large negative values of and Aq, however, one finds allowed regions in 
the parameter space with the right order of magnitude for the relic density of DM. For such 
constrained parameters, the singlino always appear to be heavy while the LSP is mainly 
bino. In some cases, the main bino annihilation channel is a pseudoscalar singlet resonance. 
We show in fig. 1(a) the various constraints in the mo, M1/2 plane for Aq = —1300 GeV 
and Ak = —1400 GeV. The large mo region is theoretically excluded as well as the region 
Ml/2 ^ 280 GeV. Along the region where the lightest stau is the LSP, at small mo, one finds 
a broad band where the DM relic density is below the WMAP upper bound. In this band, 
the bino LSP coannihilates with the stau (as in the CMSSM). Applying strictly the LEP 
limits on Higgs searches from each decay channel, as it is done in NMSSMTools, excludes a 
large fraction of the parameter space, including this WMAP compatible band. Note however 
that allowing for some theoretical uncertainty in the Higgs mass calculation (estimated to be 
~ 3 GeV) would restore most of the forbidden parameter space. In particular, most of the 
WMAP compatible bino-stau coannihilation region has rrih > 111 GeV. In the narrow region 
allowed by both theoretical and collider constraints, one finds a thin band allowed by WMAP 
where the bino LSP rapidly annihilates through a pseudoscalar singlet resonance. In fig. 1(b) 
we show the DM relic density, the pseudoscalar singlet mass and twice the bino LSP mass 
as a function of mo for the same choice of parameters as fig. 1(a) and M1/2 = 600 GeV. The 



9 



m„[GeV] (a) 

Figure 2: The WMAP allowed region in the mo, 
-1500 GeV, = -50 GeV and (b) A = .1, tan/3 = 
color code as in fig. 1 



m„[GeV] (b) 

Ml/2 plane for (a) A = .1, tan/3 = 5, = 
: 10, ^0 = -1500 GeV, = -50 GeV. Same 



sharp increase in the rehc density when 2m^ ^ mp is typical of a pseudoscalar resonance 
with a small width [23] . This is to be contrasted with the case of a narrow scalar resonance 
which shows smoother variations because the scalar does not contribute to the annihilation 
of a pair of neutralinos with zero relative velocity. 



4.1.2 Intermediate tan/3 

As explained in sec. 2, increasing tan/3 implies A < .1, larger values leading to a light scalar 
Higgs doublet excluded by LEP. We now consider A = .1, tan/3 = 5 or 10. For these 
intermediate values of tan/3, the singlino can be light if \Ao\ is close to its lower bound and 
mo, Ml/2 are small. However, for our choice of A and tan/3, small values of \Aq\ always 
lead to light states in the Higgs sector, excluded by LEP. Points for which > are also 
excluded by light Higgs states. For < 0, and large < one finds points in agreement 
with LEP constraints for which the singlino is always heavy and the LSP is mainly bino. In 
addition, small < are favorable for bino LSP annihilation through a Higgs resonance. 
The possible resonances are the pseudoscalar singlet P or the lightest scalar doublet h. As 
an example we consider Aq = —1500 GeV, A^ = —50 GeV, see fig. 2. 

The possibility of annihilation through a pseudoscalar resonance at intermediate values 
of tan/3 is a characteristic feature of the NMSSM. In our case study, the pseudoscalar singlet 
resonance is found around M1/2 ~ 300 GeV, and corresponds to |2m^ — mp] < 3 GeV. 
For larger negative values of A^, mp increases, as can be seen from eq. (2.8), which means 
that rapid annihilation would be possible for a heavier bino LSP, that is a larger M1/2. 
However, when is large it becomes increasingly difficult to rely exclusively on Higgs 
exchange to have efficient enough annihilation. Thus, for large negative values of A^, the 
rapid annihilation region disappears. The pseudoscalar singlet exchange is also dominant 
in the small region compatible with WMAP around mo ~ 2.3 TeV, see fig. 2(b). When 
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annihilation is dominated by the pseudoscalar singlet exchange, the annihilation channels 
are purely into bb and rr pairs for light pseudoscalars, whereas the tt channel can contribute 
significantly, once passed the top threshold. 

Rapid annihilation of the bino LSP through the light scalar doublet h can also occur at 
low values of M1/2 ~ 130 GeV. Note however that, as in the CMSSM, the small M1/2 region 
is constrained by chargino searches and Higgs searches at LEP, especially at small tan/3. 
For tan/3 = 5, the Higgs constraint rules out practically all the scalar doublet annihilation 
region. For tan/3 = 10, the Higgs bound is relaxed, and both the WMAP and the Higgs 
constraints are satisfied for M1/2 ~ 130 GeV. Alternatively, one could have increased \Ao\ to 
relax the Higgs bound. 

Sfermion coannihilation can also provide a mechanism to lower the relic density below 
the WMAP upper bound. As in the MSSM, we find a stau coannihilation band for values of 
mo just above the stau LSP forbidden region as well as a narrow stop coannihilation region 
at small values of mo and M1/2 just above the stop LSP forbidden region. Note that the 
latter can only be found for large negative values of Aq and it satisfies the LEP constraints 
on the Higgs sector only when tan/3 = 10. Smaller values of mo,Mi/2 are excluded as they 
lead to a negative squared mass in the stop sector. 

4.1.3 Large tan/3 

We next consider the case tan/3 = 50, A = .1. First note that one can find allowed regions for 
Aq, Aj^ < 0, as we had for intermediate values of tan/3, but also for Aq, A^ > 0. As explained 
in sec. 2, for large values of tan/3, the singlino mass is —2A\ which depends mainly on Aq but 
not on mo or Mi/2- Hence, the singlino can be the LSP for large values of M1/2, where the 
bino is heavy. However, since A = .1, small values of \Aq\ lead to light Higgs states excluded 
by LEP. Hence, \Aq\ has to be large and the regions where the singlino is the LSP are located 
at Ml/2 ~ several TeV, where no mechanism is available for singlino annihilation. 

If Aq, Ak < 0, DM annihilation mechanisms are sfermion coannihilation and Higgs ex- 
change annihilation. This is illustrated in fig. 3(a) for Aq = —1500 GeV, = —50 GeV. 
The possible Higgs resonances are again the light scalar doublet h at low M1/2 ~ 130 GeV, 
just above the chargino exclusion limit from LEP, or the pseudoscalar singlet P for slightly 
larger values of M1/2 ~ 200 GeV. For 2 < mo < 2.5 TeV and 130 < M1/2 < 200 GeV 
one also finds a WMAP allowed band where the bino LSP rapidly annihilates through the 
pseudoscalar singlet resonance. The special features of the region mo > 2.5 TeV will be dis- 
cussed in the next subsection. The LEP constraints on the Higgs sector excludes the region 
at small M1/2 where .9 < mo < 1.6 TeV. For mo ~ 750 GeV and 130 < M1/2 < 200 GeV, 
the bino LSP coannihilates with the stop NLSP. For larger values of M1/2, one also finds a 
bino-stau coannihilation thin band along the forbidden stau LSP region. Smaller values of 
mo < 750 GeV are excluded as they would lead to negative sfermion masses. 

When Aq and A^, are both positive, all Higgs states, except the scalar singlet S, can 
be light and contribute significantly to the bino LSP annihilation. In fig. 3(b), we consider 
Aq = 1500 GeV and = 250 GeV. First we focus on the region mg < 2.5 TeV. Larger values 
of mo will be treated in the next subsection. The WMAP compatible regions correspond 
to the resonances of the light scalar doublet h just above the chargino exclusion limit from 
LEP at Ml/2 ~ 130 GeV, of the heavy scalar/pseudoscalar doublets H/A in a wide band 
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Figure 3: The WMAP allowed region in the mo,Mi/2 plane for (a) A = .1, tan/3 = 50, Aq = 
-1500 GeV, = -50 GeV and (b) A = .1, tan/3 = 50, Aq = 1500 GeV, A^ = 250 GeV. Same 
color code as in fig. 1 

for 130 < Ml/2 ^ 800 GeV, and of the pseudoscalar singlet P at still larger values of 
Ml/2 ~ 900 GeV. Note however that the regions around the lighter resonances are ruled out 
by the LEP limit on the Higgs sector. The annihilation channels near the P resonance are 
into bb, TT, as well as WH^, hA or ZH. Of course coannihilation with sfermions is always 
possible when mo is near its lower bound. 

4.1.4 Large ttiq 

In the CMSSM, the interest from the point of view of compatibility with WMAP of the 
large mo region has been widely stressed. The main reason is that the parameter fi decreases 
sharply as mo increases, before entering an unphysical region where /i^ < 0. When /i < Mi, 
the LSP has an important higgsino component and can annihilate efficiently into W pairs. 

In the CNMSSM, fi also decreases at large mo, although one needs to take into account 
another factor: the squared mass of the lightest pseudoscalar Higgs state can become negative 
at large mo, thus leading to an unphysical region before the higgsino becomes LSP. This is 
precisely what happens in the example we have considered previously with tan/? = 10, in 
fig. 2(b). Large values of mo are excluded as they would lead to a negative mass squared in 
the Higgs sector and the higgsino component of the LSP remains well below 10% over the 
parameter space of the theoretically allowed region, so the annihilation into W pairs is not 
efficient and the relic density is too large. 

For tan/3 = 50 and < 0, as displayed in fig. 3(a), the higgsino component of the LSP is 
also small over the allowed parameter space at mo > 2.5 TeV. On the other hand annihilation 
through the heavy Higgs doublet H/A exchange is efficient, especially considering the tan/? 
enhanced couplings to bb and rr. This leads to an allowed band at large mo along the 
boundary of the theoretically excluded region. Note that there is a very narrow region 
at this boundary, hardly distinguishable given the scale of the figure, where the lightest 
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Figure 4: The WMAP allowed region in the mo, plane (cyan) for (a) A = .01, tan/3 = 5, 

^0 = 200 GeV, = -10 GeV and (b) A = .01, tan/3 = 10, ^0 = -20 GeV, = -50 GeV. We 
show the region excluded by theoretical constraints or by LEP searches on sparticles (black), the 
region where the stau is the LSP (blue), the region where the singlino is the LSP (grey), the LEP 
exclusion on the Higgs sector (red/hatch). In (a) we also show the LEP constraint on the Higgs 
sector for A = .001 (orange/hatch) or mt = 175 GeV (yellow/hatch) . 

pseudoscalar is excluded by LEP, just before its mass squared becomes negative. 

Finally, in the case tan/3 = 50, A^, > 0, fig. 3(b), the higgsino component of the LSP 
can be large at mo > 2.5 TeV near the theoretically excluded region, i.e. in the lower part 
of the wide WMAP allowed band. In this region, the annihilation channels are typical of a 
higgsino LSP: WW, bb as well as coannihilation channels with the charginos or with heavier 
neutralinos. In the upper part of the wide WMAP allowed band at large mo, annihilation 
is dominated by the exchange of the heavy scalar/pseudoscalar doublet H/A, as well as the 
pseudoscalar singlet P. The annihilation products include a variety of channels such as bb, 
TT, WH^, hA or ZH. Note that in fig. 3(b) it is possible to distinguish the thin region 
where the lightest pseudoscalar is excluded by LEP above the theoretically excluded region 
where its mass squared becomes negative. 

4.2 Small A: singlino LSP 

Next we focus on the regions of parameter space where a singlino LSP can be found. As 
already mentioned, this scenario appears only for A ^ 1, i.e. A = .01 in our selected scans. 
Qualitatively the results are similar for smaller values of A. In this case, we have an effective 
MSSM with an almost decoupled singlet sector and we do not expect to have singlet Higgs 
resonances. 

4.2.1 Intermediate tan/5 

While scanning over Aq, A^, mo, M1/2 with A = .01 and tan/3 = 5 or 10, we have not found 
regions where rapid annihilation through a Higgs exchange could take place. Thus the only 
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Figure 5: The WMAP allowed region in the ttiq, Mi^2 plane (cyan for singlino, green for non 
singlino LSP) for (a) A = .01, tan/? = 10, Aq = 250 GeV, = 270 GeV and (b) A = .01, tan/3 = 5, 
^0 = 750 GeV, A^ = 10 GeV (here uit = 175 GeV). Same color code as fig. 4. 

mechanism that can provide the correct relic density for a singlino LSP is coannihilation 
with a slepton NLSP. This works most efficiently when the LSP and NLSP are well below 
the TeV scale. We therefore expect to find WMAP allowed regions for choices of input 
parameters that predict a singlino LSP at low values of mo and M1/2 ^ 1 TeV. This means 
that Aq cannot be too large, cf. sec. 2. 

Let us start with < 0. As explained in sec. 2, in this case the singlino mass rrig = 2v 
is positive and grows with ttlq and Mi/2- Therefore, one expects to find a singlino LSP for 
small values of the soft masses. For tan/? = 5, = 200 GeV, A^ = —10 GeV, fig. 4(a), 
the singlino is the LSP for mo < 500 GeV, M1/2 < 1100 GeV. The singlino LSP satisfies 
the WMAP upper limit in a narrow band just below the stau LSP excluded region. There, 
the singlino LSP, which mass is rrig ~ 350 — 380 GeV, coannihilates with the stau and other 
sleptons. Most of the singlino LSP region is excluded by the LEP constraints on the Higgs. 
The whole WMAP compatible area where the LSP is a singlino can however escape the 
LEP constraints for smaller values of A (e.g. A = .001) or for a larger top quark mass (e.g. 
rrit = 175 GeV). From fig. 4(a), one can easily see the effect of A on the Higgs constraints: 
if one had taken A = .1, the whole singlino LSP region would have been excluded by LEP. 
Increasing tan/3 = 10 also reduces the impact of the LEP constraint on the Higgs mass. For 
example, assuming Aq = —20 GeV, = —50 GeV one finds a WMAP compatible singlino 
LSP region close to the the area at low mo where the stau is the LSP, see fig. 4(b). 

For Ak > 0, the singlino mass still grows with mo and M1/2 but is negative. Therefore, 
the singlino LSP is rather found at large values of mo, M1/2, just below the excluded region 
where the pseudoscalar singlet mass squared becomes negative. However, these choices 
of parameters always lead to light states in the Higgs sector, excluded by LEP. Yet for 
tan/5 = 10, Aq = 250 GeV, = 270 GeV, one still finds a region where singlino-stau 
coannihilation is just enough to satisfy the WMAP bound while LEP constraints on the 
Higgs are satisfied, see fig. 5(a). For tan/5 = 5 it is difficult to find a WMAP compatible 
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Figure 6: The WMAP allowed region in the mo, plane for (a) A = .01, tan/3 = 50, ^0 = 

-1000 GeV, A = -50 GeV and (b) A = .01, tan/? = 50, ^0 = GeV, = 50 GeV. Same color 
code as fig. 5. 



region not excluded by LEP constraints. However assuming = 175 GeV, one can find 
such an allowed region see fig. 5(b). In all cases, the mass difference between the singlino 
and the stau must be < 3 GeV for the coannihilation mechanism to be efficient enough. 

4.2.2 Large tan/3 

For tan/3 = 50, the singlino mass depends mainly on Aq, as discussed in sec. 2. Hence, for 
large values of M1/2 the bino is heavy and the singlino can be the LSP. One also expects 
more Higgs resonances as for A = .1, tan/3 = 50 above, or in the CMSSM at large tan/3. 

Let us start with the case Aq = —1000 GeV, Af^ = —50 GeV, as illustrated in fig. 6(a). 
For Ml/2 ^1-4 TeV, the situation is similar to the CMSSM: the LSP is mainly bino and its 
relic density is below the WMAP upper bound either if it coannihilates with the stau NLSP 
at low values of tuq (close to the forbidden zone where the stau is the LSP) or if it annihilates 
rapidly through a Higgs resonance. The light Higgs doublet h can play this role near M1/2 ~ 
130 GeV, although this area is mostly excluded by the LEP constraints on the Higgs sector. 
The heavy scalar/pseudoscalar doublet H/A can also play this role either for mg ~ 1 TeV, 
Ml/2 ~ 1.3 TeV or at large mo, above the excluded region where the lightest pseudoscalar 
mass squared becomes negative. When one approaches this theoretically excluded region at 
very large values of mo > 4 TeV, the higgsino component of the LSP increases and its relic 
density drops due to rapid annihilation into W pairs and coannihilation with charginos or 
heavier neutralinos. For M1/2 < 1.4 TeV, the only difference with the CMSSM is the small 
WMAP compatible zone at mg ~ 2.2 TeV, M1/2 ~ 260 GeV where the bino LSP rapidly 
annihilates though resonance of the pseudoscalar singlet P. When M1/2 > 1-4 TeV, the 
LSP is mainly singlino and its relic density is below the WMAP upper bound either if it 
coannihilates with the stau NLSP at low values of mo or if it coannihilates with the mixed 
bino-higgsino NLSP at very large values of mo > 5 TeV. Note that the WMAP compatible 
singlino LSP region at large mo does not extend all the way to the theoretically excluded 
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region: indeed, as niQ increases, /i decreases and the LSP becomes predominantly higgsino. 

For > the situation is similar. In fig. 6(b) we display our results for Aq = 0, 
Af^ = 50 GeV. Now the singlino is the LSP for M1/2 > 600 GeV and its relic density is 
below the WMAP upper bound in 3 different regions: At low mo where it coannihilates with 
the stau NLSP. For M1/2 ~ 600 GeV and 400 < mo < 1000 GeV, where it coannihilates 
with the bino NLSP which in turn annihilates rapidly through the heavy scalar/pseudoscalar 
doublet H/A resonance. Finally, at large mo the singlino LSP coannihilates with the mixed 
bino-higgsino NLSP, which in turn annihilates rapidly through the H/A resonance or into 
typical higgsino channels, including W pairs or coannihilation with charginos and heavier 
neutralinos. Here again, the WMAP compatible singlino LSP region does not extend all the 
way to the theoretically excluded region but stops where the LSP becomes mainly higgsino. 

5 Discussion 

We have achieved a first exploration of the parameter space of a constrained NMSSM from the 
point of view of DM relic density, taking into account all theoretical and collider constraints. 
We have presented our results in the mo, M1/2 plane for selected values of A, tan/5, Aq and 
A^. We have assumed sign(/i) > 0, yet we do not expect any change in our analysis for 
sign(/i) < 0. We have recovered the main scenarios of the MSSM as well as new ones. 

For A > .1 we have shown that it was possible to have an extra pseudoscalar singlet 
resonance at any values of tan/5. The search of this extra Higgs state might reveal an 
interesting challenge at LHC [67]. Although this pseudoscalar singlet state is not very 
heavy its couplings to fermions are suppressed relative to a doublet pseudoscalar making the 
pseudoscalar invisible unless tan/? is large [68]. 

For small A < .01, we have shown that it was possible to have a singlino LSP with 
a relic density below the WMAP upper bound. However, such scenarios always require 
coannihilation with the stau, bino or mixed bino-higgsino NLSP, and a small mass difference 
between the singlino LSP and the NLSP ( < 3 GeV). In these coannihilation scenarios, 
the presence of a singlino LSP in the decay of the NLSP (sfermion or neutralino) might 
influence markedly the phenomenology at colliders [54,69]. The consequences for indirect 
detection of dark matter in all the WMAP compatible scenarios will be analysed in a separate 
publication [70]. These scenarios represent only a subset of the possible scenarios for the 
singlino LSP DM in the general NMSSM. Indeed, in the model with free parameters at the 
weak scale, it is also possible to have annihilation of singlino LSP through a Z or a light 
scalar/pseudoscalar resonance or to have annihilation of a mixed singlino-higgsino LSP into 
W or Higgs pairs. Such scenarios usually require A > .1 for which the singlino is never the 
LSP in the CNMSSM^. 

We should also mention that if one assumes minimal flavour structure, the branching 
ratios for 6 — > 57 or 5^ — s> fi~^fi~ might impose strong constraints on the CNMSSM pa- 
rameter space, the latter being specially relevant for large values of tan/3 [71]. A global flt 
to all observables including those of the flavour sector in the CNMSSM is left for a future 
publication. 

■^Previous phcnomenological studies of the NMSSM with strict universahty at the GUT scale also con- 
cluded that A < .01 was required in order to have a singlino LSP [54]. 
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